Functional integral approach for multiplicative stochastic processes 
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We present a functional formalism to derive a generating functional for correlation functions of 
a multiplicative stochastic process represented by a Langevin equation. We deduce a path integral 
over a set of fermionic and bosonic variables without performing any time discretization. The usual 
prescriptions to define the Wiener integral appear in our formalism in the definition of Green func- 
tions in the Grassman sector of the theory. We also study non-perturbative constraints imposed by 
BRS symmetry and supersymmetry on correlation functions. We show that the specific prescription 
to define the stochastic process is wholly contained in tadpole diagrams. Therefore, in a supersym- 
metric theory the stochastic process is uniquely defined since tadpole contributions cancels at all 
order of perturbation theory. 
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I. INTRODUCTION 

Stochastic processes [l[ has had a great interest in the 
scientific community since many years ago. Applications 
can be found along a wide area of scientific research, from 
physics and chemistry 0, Q through biology and ecolo- 
gy HLHl to economy and social sciences @,I3- 

Perhaps, the most popular formalisms to deal with 
this type of problems are the Langevin and Fokker-Plank 
ones 0, [|[ . While the Langevin description is based on a 
set of stochastic differential equations, the Fokker-Plank 
one is based on the probability distribution of the set 
of variables satisfying a deterministic partial differential 
equation. The Langevin description rose as a formula- 
tion of Brownian motion [8[ and gives a very clear in- 
sight into the microscopic properties of the system. In 
its original application [9j, the influence of the medium 
on the diffusive particle is modeled by splitting its effect 
in a deterministic part, given by the viscous force, and 
a stochastic part, given by a random force with zero ex- 
pectation value. In this way, fluctuations exhibit as an 
additive noise and, consequently, the considered model is 
an additive stochastic process. In a different case, fluctu- 
ations could depend on the state of the system [l(| and 
they are regarded as the product of the random force and 
a function of the state variable. In that situation, mul- 
tiplicative noise is defined and the model is known as a 
multiplicative stochastic process. 

There is a continuously growing interest in modeling 
and studying noisy systems j7l| and a great amount of 
work was developed for additive as well as multiplica- 
tive noise. As an example of multiplicative noise, we 
would like to mention the stochastic Landau-Lifshitz- 
Gilbert [l2j equation, used to describe dynamics of classi- 
cal magnetic moments of individual magnetic nanoparti- 
cles. Here, the noisy fluctuations of the magnetic field 
couple the magnetic moment in a multiplicative way. 
Another interesting example is the problem of state- 



dependent diffusion [lfjj. The simplest example is the 
diffusion of a particle in suspension near a wall where 
the diffusion coefficient depends on the distance to the 
wall. 

We are interested in the path integral formulation of 
stochastic processes modeled by Langevin equations with 
multiplicative noise. This is an interesting topic as long 
as the path integral formalism provides a useful tech- 
nique to compute correlation and response functions [l3| . 
Once one is able to built a generating functional for cor- 
relation functions, it is possible to study the system us- 
ing all the machinery developed in quantum field the- 
ory, even non-perturbative techniques. In particular, it 
would be possible to implement a dynamical renormal- 
ization group [id], EH to look for the asymptotic behavior 
of correlation functions. Within this formalism, symme- 
tries play an essential role since they give information 
about non-trivial conserved quantities. Moreover, non- 
trivial symmetries [16| govern the route to equilibrium of 
a dynamical system. 

The first quantum field theory formalism describing 
additive noise was developed by Martin, Siggia and 
Rose [l7j to perturbatively compute correlation and re- 
sponse functions. After that, the formalism was exten- 
sively developed and today, the case of additive noise 
is reasonably well understood [Hj]. Path integral meth- 
ods for multiplicative noise have been derived in different 
situations. When the stochastic process is driven by col- 
ored noise there is no interpretation problem and a phase- 
space path integral treatment was dev elop ed for Gaussian 
as well as non-Gaussian shot noise [19j . For Gaussian 
white noise, path integral methods were worked out in 
the Ito [13], Stratonovich [2l| and also in more general 
conventions [IcJEIj]- However, the case of multiplicative 
noise is more involved and some controversies [22[ have 
appeared in the literature. Indeed, for fixed stochastic 
prescription, different results were presented for the gen- 
erating functional [13, EE EE [HJ , depending on the use 
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of continuous or discrete time methods to derive the path 
integral. 

In this paper, we present a derivation of the gener- 
ating functional using just functional methods, without 
discretizing the Langevin equation at any intermediate 
step. We treat additive as well as multiplicative noise 
in the same footing. The use of auxiliary Grassman vari- 
ables allows us to describe (at least formally) a stochastic 
multiplicative process without any reference to the def- 
inition of the Wiener integral. In this way, the known 
prescriptions [23[ to define the stochastic process, emerge 
in our formalism in the definition of Green functions at 
definite points. This procedure is useful to study time re- 
versal symmetry and, in general, to study the approach 
to equilibrium. 

In the case of additive processes, there is no ambigui- 
ties in the definition of the Wiener integral. However, the 
explicit construction of a generating functional does de- 
pend on the type of discretization adopted [24] . We show 
that there is no contradiction since, due to a hidden sym- 
metry (supersymmetry) , the tadpole contributions (car- 
rying the discretization ambiguities) cancel to all order 
of perturbation theory. In addition, supersymmetry en- 
force time reversal symmetry and fluctuation dissipation 
relations. This observation may be useful to analyze non 
trivial multiplicative vector processes. We discuss in de- 
tail the relation between supersymmetry, equilibrium and 
the Generalized Stratonovich convention [23[ (also called 
"a-convention" [l5j in the field theory literature) in ad- 
ditive as well as multiplicative stochastic processes. 

The paper is organized as follows: In section [TT1 we de- 
fine our model by means of a Langevin equation with 
multiplicative noise and derive the generating functional 
for correlation functions. This is the main result of the 
paper. In ^Illl we show, carefully integrating over Grass- 
man variables, the connection of our result with previous 
ones obtained by discretization of the Langevin equation. 
Then, we study the symmetry properties of our action in 
section HVl In particular, we show that in the case of ad- 
ditive noise, where the correlation functions are indepen- 
dent of particular prescriptions, tadpole diagrams cancel 
to all order of perturbation theory due to a hidden su- 
persymmetry. Finally, we discuss our results and future 
perspectives in section [V] Appendix SfX] briefly summa- 
rizes the usual operations with Grassman variables used 
along the paper. 



II. GENERATING FUNCTIONAL 

For simplicity, we consider a single random variable 
x(t) satisfying the first order differential equation 



dx(t) 
dt 



= f(x(t)) + g(x{t))r){t), 



(1) 



where rj(t) is a Gaussian white noise, 

(»?(*)> = , ( V (t)ri(t')) = S(t-t'). (2) 



The generalization to a system of variables is straightfor- 
ward. In Eq. (TT]), the drift force f(x) and the square root 
of the diffusion function g{x) are, in principle, arbitrary 
smooth functions of x(t). 

To completely define eq. ([T]), it is necessary to give 
sense to the ill-defined product g(x(t))rj(t), since r](t) 
is delta correlated. The underline microscopic physics 
should fix this interpretation using the Ito, Stratonovich 
or even more general prescriptions. Since our formalism 
is not based on an explicit discretization of the Langevin 
equation, we will formally work independently of any 
kind of prescription. In fact, the main result, eq. (|18p . 
does not depend on the explicit definition of the Wiener 
integral. Of course, the ambiguities will show up in a 
different way, since they are at the stem of the stochas- 
tic process. We discuss this problem with great detail in 

am 

We are interested in computing n — points correlation 
functions (x(ti) . . . x(t n )). To compute them, we should 
know the n th -order joint probability of the random vari- 
able x(t). An alternative equivalent procedure is to solve 
the Langevin equation (Eq. ([1}) and compute, 



{x(ti) . ..x(t n )) = (xwih) . ..x [v] {t n )) j 



(3) 



where x\ v ](t) is a solution of Eq. (fTJ) for a particu- 
lar realization of the noise and certain initial condition 
x(to) = xq. < . . . >rj means the stochastic mean value in 
the rj variable. If the initial condition is not determine 
but it is distributed with certain probability V(xq), the 
expectation value should be taken also in this variable. 

Correlation functions can be obtained from a generat- 
ing functional 



Z[J(t)] = / e /* J W*wW 



(4) 



by simply differentiating with respect to the source J(t), 



5 n Z[J] 



5J(t n )...6J(h) 



(5) 



.7=0 



(note that Z(Q) = 1 by definition). 

The aim of this section is to find a functional represen- 
tation for Z( J) avoiding the problem of explicitly solving 
the Langevin equation. For this purpose we introduce 
a functional integral over x(t) and a delta-functional, 
rewriting Eq. (j4|) in the following form, 

Z[J{t)\ = Q Dx(t) S[x(t) - x [v] (t)} e /*'(*M*^ (6) 

or, since the only noise dependence comes from X[ v ](t), 

Z[J(t)} = ( Dx(t)eS dtJ W*W(6[x(t)-x [ri] (t)}) (7) 



The following property of the delta-functional is very 
useful, 



S[x(t) - x [v] (t)} = 6[0(x)] det 



50_ 

8x 



(8) 
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where xt v ](t) is a root of O(x), i.e., 0(x\ n ]n\) = 0. 
From Eq. ([IJ we can choose for 0{x) 



0(x) 



dx(t) 
dt 



f(x(t)) - g(x(t)) V (t) 



(9) 



and correspondingly, the differential operator SO/Sx is 
given by 



SO{x(t)) 
8x(t') 



dt 



f(x) - g'(x)r,(t) 



6(t-t'), (10) 




where /' and g' mean differentiation with respect to x. 
Therefore, Eq. ([7]) now reads 



with the definitions Eqs. © and ([TO]). 

It should be noted that for additive noise processes, 
g ' (x) = and then, the determinant in Eq. (|11[) is noise 
independent. However, in the general case, the deter- 
minant should be carefully treated. To do that, we in- 
troduce a couple of Grassman functions and 
obeying the usual Grassman algebra 

= = = o, (12) 

where { , } represents anticommutators. These relations 
imply, in particular, £(i) 2 = £(i) 2 = 0. The determinant 
can be represented [ISJ as (for details see Appendix fAf : 



det 



SO 



= j D^D(, e i dt ^ {t) ^ m e ^-l' dt ^ t ' > \- } " {x)+9 ' (x)r,{t ^ m . 

(13) 

Note that, different from the additive case, the noise cou- 
ples to the Grassman variables. 

A last auxiliary function ip(t) is introduced to represent 
the delta-functional 

I 



Substituting the results (fTSJ and (H4J) into Eq. (0, we 
find 



(8[x(t)-x [n] (t)]\ = J D£D£Dcp x 

x e fdt {mat)-i(tK(t)f^)-Mt)[i(t)-f(x)]} 

- f dt {g'(x)£{t)£(t)-iip(t)g{x)}ri(t) 



x { e 



(15) 



where the "dot" means time differentiation. 

With this representation for the delta functional, the 
average over the noise can now be easily computed. We 
find, 

(^ e - 1 dt {g K-ivg}n^ _ e ~ f dt {hf 2 9 2 +i<pgg'ii\ _ £jg) 

The first term, proportional to g 2 {x), is the usual term 
which comes from the average of the delta-functional. It 
is also obtained when dealing with an additive process. 
However, the second term, proportional to g(x)g'(x), is 
proper of multiplicative noise processes and has its ori- 
gin in the noise dependence of the determinant, stated 
in Eq. (|5|). This term is responsible for the well-known 
spurious drift terms[23]. One important fact in the cur- 
rent example of a single random variable is the absence 
of any term proportional to g'{x) 2 . It is automatically 
canceled due to the property of the Grassman variables, 
£,(t) 2 = £(£) 2 = 0. This situation changes for a system of 
stochastic equations with several variables. 

Introducing Eq. (jT5j) into Eq. (JTTJ) and using Eq. (TTo) 
we finally arrive at the desire representation for the gen- 
erating functional of correlation functions, 



Z[J] = J DxDipDiDi e S[x,<p,a]+f dtJ(t)x(t) 



(17) 



where the "action" S is given by 



dt 



-m^ffl + fwmm + \^(t) 2 g(,x) 2 + up® 



dx 



■g{x)g\x)mm 



r 



(18) 



Thus, we have deduced a generating functional for cor- 
relation functions as a functional integral over the vari- 
ables x(t), ip(t),£(t),£(t). The multiplicative noise effects 
are evident in the coupling of the Grassman variables 
with <p(t). This coupling is absent in the additive case. 



Usually, the effect of the determinant in Eq. (JSj) have been 
treated by discretizing the Langevin equation. In this 
context, there are some ambiguities in taking the contin- 
uum limit that provide, in general, different stochastic 
evolutions. In the next section, we will show how these 
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ambiguities appear in our formalism and the relation of 
our results with previous ones. 



III. INTEGRATION OVER GRASSMAN 
VARIABLES, AMBIGUITIES AND 
PRESCRIPTIONS 

We have obtained the generating functional and its 
corresponding action, Eqs. (fl7|) and (|18[) . by means of 
functional methods, without any reference to discretiza- 
tion of the Langevin Eq. (JXJ) . However, it is well known 
that different forms of discretization and therefore, dif- 
ferent ways to derive the continuum limit, produce differ- 
ent stochastic evolutions. The reason for that behavior 
resides in the delta-correlated noise. In fact, for white 
noise, the product g(x(t))i](t) (in Eq. ([1])) is ill-defined 
due to the infinite variance of the noise. In many physical 
applications this problem can be overcame by consider- 
ing a weakly colored Gaussian-Markov noise with a finite 
variance [l9j. In this case, there is no problem with the in- 
terpretation of eq. (p} and we can take the limit of infinite 
variance at the end of the calculations. This regulariza- 
tion procedure is equivalent to the so called Stratonovich 
interpretation 0, fl8| . However, in other applications, 
like chemical Langevin equations or econometric prob- 
lems [alii], the noise can be considered principally white, 
since it could be a reduction of jump-like or Poisson like 
processes. In such cases, the ltd interpretation should 
be more suitable. Therefore, the interpretation of eq. ([T]) 
depends on the physics behind a particular application. 
Once the interpretation is fixed, the stochastic dynamic 
is unambiguously defined. 

Our formalism is not obviously depending on such in- 
terpretations (see eq. (fT8|l). In this way, it is specially 
appropriated to study dynamical properties that are in- 
dependent of particular prescriptions. In the following, 
we show how different interpretation of eq. ([TJ appear in 
the Grassman path integral formalism. 

The problem can be easily understood looking at the 
integral 



the prescription for the choice of Tj. There are several 
prescriptions to define this integral that can be summa- 
rized in the so called Generalized Stratonovich prescrip- 
tion [23[ or "a-prescription" for which we choose 
g((l — a)x(tj) + ax(t,j + i)) with < a < 1. In this 
way, a = corresponds with the Ito interpretation and 
a = 1/2 coincides with the Stratonovich one. Moreover, 
the post-point prescription a = 1 is also known as the ki- 
netic or the Hanggi-Klimontovich interpretation [23|, [25l - 

In this section, we show that the same ambiguities ap- 
pear in the action, Eq. (|18p. without any reference to 
discretization. The problem is associated with the defi- 
nition of the determinant of Eq. ([T5]) . To compute this 
formal expression is necessary to define the Green func- 
tion of the operator d /dt. This problem is common to 
any integration over Grassman variables and it is not 
proper of multiplicative noise. As we will see it is also 
present in any white noise stochastic process. 

For concreteness and to illustrate this point, let us con- 
sider the following functional integration over Grassman 
functions, 



1 G = / D£,DH e 



/dtf(t)A£(t) - / dt £(t)«(t)Mt) 



(21) 



where A is any differential operator and h(t) is an ar- 
bitrary function of time. To compute it, we perform a 
Taylor expansion of the second exponential, obtaining 



x ji-y dt£(t)t(t) h(t)+ 

+ ~ J dtdf mmat'W) h(t)h(t>) - . . J .(22) 

Computing the Gaussian integral by using Wick theorem 
we find (see appendix 



In = Det 



g(x(t)) v (t)dt = I g(x(t)) dW(t) , (19) + y / dtdf [(A^iO)) 2 - A~ 1 (t,t')A~ 1 (t' ,t)] h{t)h(t') 



where we have defined the Wiener process W(t) as r)(t) 
dW{t)/dt. 

By definition, the Riemann-Stieltjes integral is 



dtg(x(t))dW(t) = lim Y,9(^))(W(t 3+1 yW(t 3 )) 

i=i 

(20) 

where Tj is taken in the interval [tj,tj + {\. For a smooth 
measure W(t), the limit converges to a unique value re- 
gardless the value of Tj. However, W(t) is not smooth, 
in fact, it is nowhere integrable. In any interval, white 
noise fluctuates an infinite number of times with infinite 
variance. Therefore, the value of the integral depends on 



- ^- / dtdt'dt" 
3! 



(23) 



where the Green function A 1 {t 1 t l ) is defined as 
A A _1 (i,t') = S(t - t'). 

At this point it is necessary to precisely define the 
Green function A _1 (t — t'). We choose a causal prescrip- 
tion by considering the retarded Green function in such 
a way that 



A- R \t-t')A- R \t'-t) = Q 



(24) 



except, possibly, in the null measure set t = t' . 

To go forward, it is necessary to carefully analyze the 
function h(t). For delta autocorrelated functions the 
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second line of Eq. (123)) and all the subsequent terms 
are automatically zero. This is a natural consequence 
of the anticommuting character of the Grassman vari- 
ables. In this case, the result is linear in h(t). However, 
if h(t) is a smooth function of time we can safely put 
A' x (t - t^A- 1 ^' - 1) = in Eq. (J23j) since t = t' is a 
null measure set. Then, it is possible to re-exponentiate 
the Taylor expansion, obtaining 



T G = Det 



-) ex P {A 



fl(0) 



dt h(t)}. 



(25) 



So, the result of Eq. (|2"5)l is valid provided we choose 
a causal prescription for the Green function of the op- 
erator A and we consider smooth functions h(t). Even 
in this situation, there is an ambiguity, given by A^(0), 
which has not been defined. To completely define the 
determinant we need to fix this value. 

Let us apply this result to integrate the Grassman vari- 
ables in the generating functional, Eq. (IT7| . In this case, 
the operator A = 4r. Therefore, the retarded Green 
function is A^ 1 — 9(t — t'), where 9(t) is the Heaviside 
distribution. Therefore, to define the integration we need 
to fix 9(0) = a. With this choice, the generating func- 
tional reads, 



Z[J] = J DxDip e s[xM+S dtJ{t)x{t) ; ( 26 ) 
where the "action" S[x,ip] is given by 



S = I dt ^ —<p g + iip 



af\. (27) 



dx ; 
-dt~ f + a " 



Equation (|27|) coincides with that of ref. |10( , where the 
parameter a was introduced discretizing the differential 
equation. Also, our action (with a = 1/2) coincides with 
the one presented in ref. |22j |. However, it differs from 
that of refs. [TEGl. 

To conclude this section, we note that the informa- 
tion about the particular interpretation of eq. (fTJ) is not 
contained in the action, eq. (|T£1) . It is encoded in the 
definition of the Grassman Green functions at the ori- 
gin. This is a general result and does not depend on the 
details of the stochastic differential equation. 

We would also like to remark that , even in the case 
of additive noise (</ = 0), the action does depend on a 
through the last term in eq. (f2"Tl) . However, Z(J) should 
be a independent (dZ/da = 0). Therefore, it should 
be a non-trivial symmetry of the action enforcing this 
property. We address and detail this issue in the next 
section. 



IV. SYMMETRIES AND CONSTRAINTS 

The introduction of Grassman variables allows to study 
symmetries, implemented as linear transformations be- 



tween the variables x, <p, £, £. These symmetries have di- 
rect physical consequences and they strongly constraint 
the correlation functions, giving information about non- 
perturbative structure of the theory. Even in the case of 
perturbation theory, they are useful by helping to rear- 
range diagrams more efficiently. 



A. BRS symmetry 

By construction, there is a hidden symmetry in the 
generating functional, Eq. (fTTj) . that is related with the 
fact that Z(0) = J Vr] V(r]) = 1 (see Eq. g))), i. e., it is 
a consequence of probability conservation. 

From Eq. (jlip it can be readily seen that 

Z[0] = l^J Dx(t) det (^j 5[0(x)]j 

= Q DdS[6}^ =(!), = !, (28) 



where the measure 



DO = Dx(t) det ( ^ 



(29) 



is trivially invariant under the translation group O — > 
O + a, with constant a. Since we are working with the 
stochastic variable x(t) 1 this symmetry is implemented as 
0(x') = 0(x) + a, that, for infinitesimal transformation, 
reads 



Sx = x — x = 



SO 
Sx 



(30) 



where a is an infinitesimal parameter and ^^rl is the 

inverse operator of Eq. (fTU| . Thus, the invariance of the 
measure, Eq. (|29|). under the highly non-linear transfor- 
mation Sx, guarantees the probability conservation prop- 
erty Z(0) = 1. In this way, it is clear the role played by 
the determinant as a Jacobian of the non-linear change 
of variables O = 0(x). 

The representation of the determinant in terms of 
Grassman variables allows to visualize the non-linear 
transformation, Eq. (|30l) . as a linear transformation be- 
tween variables in the extended space (x, cp,£, £). Indeed, 
it can be shown that, even after the integration over the 
noise, the generating functional, Eq. (flTl) . is invariant 
under the linear transformations 



Sx = A£ , S£ = , 
<5£ = i\ip , Sip = , 



(31) 
(32) 



where A is an anticommuting parameter. This nilpotent 
transformation (S 2 — 0) is the famous BRS [28| symme- 
try, discovered in the context of quantization of Gauge 
theories. 
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A more convenient representation is obtained by intro- 
ducing a new Grassman ("temporal") coordinate and 
the functions 

X(t,9) = x(t)+9£(t), (33) 
C(t,9) = t(t)+i9<p(t). (34) 

In terms of these functions, the transformations in Eqs. 
(f3~Tj) and (j3"2")l are simply translations in 0. In fact, 

SX(9) = A VX{9) , SC(9) = A VC{9) , (35) 

where the generator V = 8/ 89, satisfying T> 2 = 0. 

The action of Eq. (IT51) can be written in terms of X and 
C in such a way that the BRS symmetry displays in an 
obvious way. We can write the action, before integrating 
over the noise, in the following way: 

S = J dtd9 0(6) (x(9) - f(X(9)) - g(X(9)) r,(t)) . 

.(36) 

It is interesting to note that, in the generating functional, 
the formal functional integration over C enforce the equa- 
tion 

X{9) - f{X{9)) - g{X{9)) V (t) = 0, (37) 

which has the same (formal) structure of the Langevin 
equation, Eq. ([1]). 

Integrating over the noise, we find 

S = J dtd9 (7(0) (x(9) - f(X(6)f) 

+ J dtdhdhc{h)c{e 2 )g{x{e x ))9{x{e 2 ))-{^) 

The advantage of this expression is that it is explicitly 
BRS invariant, since its dependence on 9 is only implicit 
thorough the functions X and C. Of course, integrating 
Eq. J3HD over d9, we find Eq. (JTBJ) as we should. 

From the representation of the determinant, Eq. (IT3"1) . 
we observe that the Grassman variables always appear in 
pairs ££. Assigning a "fermionic number" +1 to the vari- 
able £ and — 1 to £, each term of any correlation function 
must have zero fermionic number, i. e., there is fermionic 
number conservation. This can be seen in eq. (|38[) by the 
fact that for each d9 there is one function C. Fermionic 
number conservation imposes 

Wx)x{t2)) = (€(*i)y»(*a)> - (C(*i)e(*2)> = 0, (39) 

therefore, (X(t x , 9~i)X(t 2 , h)) = {x(t{)xfo)). 

On the other hand, since the system is invariant under 
translations in 9, any correlation function should be a 
function of 9 2 — 6\. For instance, the correlation (XC) 
takes the general form 

(X(ti,fli)C(t 2 ,fl 2 )) = A(t u t 2 ) + (§ 2 - X ) B{h,t 2 ), 

(40) 



where A and B are functions of time. Combining this 
equation with the definitions of X and C, Eqs. (|33j) - 
(f34|) , and considering the fermionic number conservation, 
Eq. ([23), we find 

(x(hMt 2 )) = (^i)f(ia)), (41) 

in such a way that 

(x(h, 9x)C(t2, e 2 )) = (e 2 - h) (x(txMt 2 )). (42) 

Eq. (|41j) is an important constraint imposed by BRS sym- 
metry and it is valid to all order of perturbation theory. 
It states that the correlation functions of Grassman vari- 
ables are equal to the physical response of the system to 
a delta perturbation. 

We have seen in the last section that, in order to in- 
tegrate over Grassman variables, it is necessary to de- 
fine the Green function. We chose the retarded prescrip- 
tion, supplemented with A~j^(t,t) = a. Eq. (j4T|) imposes 
the same behavior for the response functions, in fact, 
(x(tx)tp{t 2 )) = for tx < t 2 and (x(t)<p(t)) = a. Thus, 
we are not free to choose any other prescription for the 
responses. If, for some reason, we decide to use another 
convention we would be breaking BRS symmetry and, in 
that way, we would be spoiling the probability conserva- 
tion. 



B. Supersymmetry and tadpole cancellation 

There is another convenient representation of the same 
action that reveals new symmetries and constraints. Let 
us introduce a new Grassman variable 9, conjugate to the 
previous one 9 and define the superfield 

$(t,9,§) = X(t,§) + C(t,§)9. (43) 

It is simple to see that 

c(9) = -|| , f(x(6)) = J deem. (44) 

Then, Eq. (|3"8)) can be rewritten as a local action in the 
superspace {t, 9, 9} 

J [ 89 \ 8t V ' } 2 89 d9 y V ' 

(45) 

In this representation, BRS symmetry is again evident 
since the action does not depend explicitly on 9 (note 
the asymmetric role of 9 and 9). The correlation func- 
tion 9x, ti)Q(9 2 , 9 2 , t 2 )), which contains the phys- 
ical correlation functions between x(t),<p(t),£(t),£(t), is 
strongly constrained. Due to the BRS symmetry, causal- 
ity and fermionic number conservation, a general two 
point correlation function, has the form 

($i$ 2 ) = {x(t 1 )x(t 2 )} + (9 1 -9 2 )x 
x [(0i - 9 2 ){x(t 1 )i P (t 2 ))- + (0i + 9 2 )(x(txMt 2 )}+] , 

(46) 
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where =F represents retarded and advanced prescriptions, 
respectively. 

Taking into account that the parameter a = 
(x(t)tp(t)) = (£,(t)(,(t)), any correlation function depends 
on a only through tadpole diagrams in a perturbativc 
expansion. In the superspace formulation, tadpole dia- 
grams represent products of superfields and their deriva- 
tives at the same superspace point (t, 9, 9). 

Interestingly, from Eq. (l46l) it is immediate to check 
the property [29| 



dtdddd 9, 0)$(t, 6, 9)) = , 



(47) 



which has deep consequences in the dynamics of the 
system. To see this point, let us consider the exam- 
ple of a Langevin equation with additive noise, i.e., 
<7($) = g = constant. In this situation, there is an- 
other BRS symmetry related with the variable 9, which 
constraints even more the correlation functions. In fact, 
it is possible to rewrite the action, Eq. (|45l) , to make this 
symmetry explicit. In the case of one random variable, 
it is always possible to write /($) = dV Jffl ■ Then, 



S = J dtd9d9 + U($)} , 



where we have defined 



D 



d_ 

00 



1 ,<9 
D =2 9 09 



dt 



(48) 



(49) 



which satisfy D 2 = D 2 = and {D, D) - 



dt ■ 



The action of Eq. (|48p has an extended symmetry 
(called supersymmetry (SUSY) 30]), whose generators 



are 



D' 



d_ 
86 



— . 1 ? d 

= 2«m 



dt 



(50) 

The first generator, D' coincides with the previous BRS 
one. In addition, there are two more generators asso- 
ciated with translations in time and in 9. The graded 
algebra 

{D',D} = {D',D} = {D',D} = {D',D} =0,(51) 
D' 2 = D' 2 = , (52) 

guarantees that the variation of the action is a total 
derivative. Supersymmetry has very important physical 
as well as mathematical implications. From the phys- 
ical point of view, it ensures that the stochastic pro- 
cess reaches the thermodynamical equilibrium at long 
times [l6l. Il8l|. In fact, SUSY implies that two-point cor- 
relation functions depend on the time difference. On the 
other hand, it enforces the fluctuation dissipation theo- 
rem [13. 

From the mathematical point of view, noting that 
V(3>) does not depend on "velocities" (nor time, neither 
9 derivatives) and using the general equation Eq. (|47[) . it 



is immediate to conclude that the contribution of tadpole 
super-diagrams cancels identically to all order of pertur- 
bation theory. This is a general result in any supersym- 
metric theory [3(| • It is based on the fact that the Hamil- 
tonian (i. e., the generator of time translations, H) can 
be cast as an anticommutator of SUSY generators and, 
since the ground state is supersymmetric, (H) = 0. 

As long as there are no more diagrams, other than tad- 
pole ones, which contain (x(t)ip(t)) — (£(t)£(i)) = a, we 
conclude that correlations do not depend on the prescrip- 
tion to define Green functions of differential operators. 
SUSY guarantees that the stochastic dynamic is unique 
and tends to equilibrium configurations at long times, 
independently of any prescription to define the Wiener 
integral. Therefore, correlation functions are a indepen- 
dent at any order of perturbation theory. 

In the multiplicative noise case, we have no such a 
symmetry for general functions f(x) and g(x) (Eq. (|45p). 
Although Eq. (|4T|) is valid, because it only depends on 
BRS symmetry, there will appear tadpole diagrams with 
derivative couplings like ) that do not cancel and 

therefore, correlation functions will depend on a. In this 
situation, each prescription a defines a different stochas- 
tic process. Thus, the coupling of the Grassman variables 
with the stochastic variable x(t) through gg'££ gives rise 
to the usual spurious drift terms. It should be noted 
that, in the case of several variables (vector variables), 
the spurious drift term could be zero for a special class 
of diffusion matrices (2^. In this case, the correlation 
functions do not depend on a even for the case of mul- 
tiplicative noise. Whether this behavior implies a super- 
symmetry in the path integral approach is an interesting 
open question. 



V. DISCUSSION AND CONCLUSIONS 

We have developed a functional formalism to describe 
stochastic processes modeled by Langevin equations with 
multiplicative noise. We have treated additive as well 
as multiplicative noise at the same level. Our deriva- 
tion of the generating functional for correlation and re- 
sponse functions does not rely on the discretization of 
the Langevin equation at any intermediate step. In this 
way, we did not need to explicitly face the ambiguities in- 
herent to the definition of a Wiener integral. Of course, 
since these ambiguities are at the stem of the stochas- 
tic process, they reappear in our formalism in a different 
context, i.e., in the precise definition of Green functions 
in the Grassman sector. 

We found a generating functional written in terms of 
a functional integral over bosonic as well as fermionic 
(Grassman) variables (Eqs. (fT?) and (fISj)). We have 
also shown that, carefully integrating over the Grassman 
variables, we reproduce previous results [Io|, HH obtained 
by discretizing the Langevin equation (Eq. (|27p). 

In particular, we have shown that the results achieved 
by using the "a-convention" are reobtained in our for- 
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malism by just fixing the retarded Green function for 
Grassman variables, Gji(t,t) — a. In our example, the 
differential operator is 4 and the retarded Green func- 
tion is the Heaviside distribution 0(t — t'). However, the 
method is general and can be applied to any differential 
operator. 

One of the practical advantages of this formalism is 
that it is not necessary to take care of the effects that 
prescriptions make on different stochastic evolution. It is 
possible to work with the same mathematical structure 
and, at the end of the calculation, the equal time Green 
function is fixed, whether it is necessary. 

Also, the presence of Grassman variables makes the 
symmetries of the theory evident. When the action is 
formulated in superspace, in terms of 0, 0), it is clear 
that the stochastic prescriptions are related with tadpole 
super-diagrams. Moreover, if the theory is supersymmet- 
ric, the tadpole diagrams cancel identically to all order 
of perturbation theory. Therefore, in the special case of 
a supersymmetric theory, all prescriptions to define the 
Wiener integral conduce to the same stochastic evolution. 
On the other hand, from a physical point of view, su- 
persymmetry implies that the stochastic evolution tends 
asymptotically to the thermodynamic equilibrium [32j. 

We believe that our approach clarifies previous con- 
troversies [22[ about the difference between discrete and 
continuum approaches to the path integral formulation 
of multiplicative noise. On the other hand, it provides a 
natural unified framework to deal with different prescrip- 
tions used to define stochastic processes. The implemen- 
tation of BRS symmetry and supersymmetry in a linear 
form allows to compute restrictions to correlation func- 
tions in a non-perturbative way, helping to guide actual 
computations in concrete models. 

The formalism presented in this paper can be readily 
generalized to more general cases as, for instance, to vec- 
tor stochastic variables. It could also be possible to treat 
systems driven by white shot noise. We hope to report 
on those issues briefly. 
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Appendix A: Representation of a determinant in 
terms of Grassman integration 

We briefly review in this appendix some usual manipu- 
lations with Grassman variables that we have used along 



the paper. For a detailed treatment of this subject in 
statistical mechanics as well as in quantum field theory 
we refer the reader to ref. 

We define a set of n Grassman variables 0i and its 
conjugates 0j, where i — 1 . . . n as 



{9i,6 j } = {6i,§ j } = {0 i ,9 J } = O 



(Al) 



where {a, b} = ab+ ba is the anticommutator of a and b. 
This definition implies the nilpotent property 0\ = 0\ = 
0. Therefore, any function of these variables should be a 
polynomial of degree at least In. For instance, if n = 1, 



F(0, 



A+B0+C9+D 



(A2) 



is the most general function of and 0. The coefficients 
A, B, C, D are complex numbers. 

Differentiation in the Grassman variable is also a nilpo- 
tent operator d 2 /d 2 8 — satisfying the Clifford algebra 



{ 



d_ d_ 

del' M 



} = o, { 



d_ 

eel 



(A3) 



and its conjugates. Any anticommutator that mixes 9i 
and 9j is also zero since they are independent. 

Interestingly, integration in a Grassman space is the 
same operation as differentiation. Therefore, the use of a 
derivative or an integral symbol is a matter of taste. For 
instance, taking into account eq. (|A2[) . 



J d0d0 F(0, 0) = D (A4) 

That means that the integration over d9d9 picks up the 
coefficient of the term 00. 

Let us consider the Gaussian Grassman integral 



Ig{A) 




(A5) 



fc=i 



where j4y are the elements of a matrix A. According 
with eq. (IA4|) . the integral is the coefficient of the term 
proportional to 0\0\0i0i ■ ■ ■ 9 n 9 n - Therefore, expanding 
the exponential in eq. (| A5|) in a "finite" Taylor series, 
and reordering the terms taking into account the anti- 
commuting properties of the Grassman variables, 

I G (A)= (-l) P A ndn A n _ ldn _ 1 ...A ldl (A6) 



where P is the order of the permutation of {j-y, ... ,j n }. 
We immediately recognize that 



I G (A) = det(A) 



(A7) 



This result should be compare with the output of a nor- 
mal Gaussian integral that is 1(A) ~ det^ 1 (A). There- 
fore, in the same way that the inverse of an n x n matrix 
determinant can be represented as a Gaussian integral 
over a set of n complex variables, the determinant itself 
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can be represented as a Gaussian integral over a set of n 
complex Grassman variables. 

We can generalize now the case of a discrete set of 
Grassman variables {6*i, . . . 9 n }, to an infinite set of con- 
tinuous variables; i. e. , a Grassman function £(t). In 
this case we can generalize eqs. (|A5|) and (|A7|) to 

det(A)= [ DtDgef (A8) 



where det(A) is a functional determinant and D^D^ are 
functional integrations over Grassman variables. A(t,t') 
is the kernel of the functional A. This is the formula used 
in eq. (| 13[) to represent the Jacobian 50/5x. 
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